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The evolutionary significance of stress-induced mutagenesis was evaluated by 
studying mutagenesis in aging colonies (MAC) of Escherichia coli natural iso- 
lates. A large fraction of isolates exhibited a strong MAC, and the high MAC 
variability reflected the diversity of selective pressures in ecological niches. 
MAC depends on starvation, oxygen, and RpoS and adenosine 3’,5’-monophos- 
phate regulons; thus it may be a by-product of genetic strategies for improving 
survival under stress. MAC could also be selected through beneficial mutations 
that it generates, as shown by computer modeling and the patterns of stress- 
inducible and constitutive mutagenesis. We suggest that irrespective of the 
causes of their emergence, stress-induced mutations participate in adaptive 


evolution. 


Bacteria are champions of evolutionary suc- 
cess; they grow in practically all ecological 
niches. Evolutionary success depends on phe- 
notypic selection, which in turn depends on 
available genetic variability. The genetic 
variability is produced primarily by mutagen- 
esis and secondarily by recombination, which 
shuffles preexisting mutations. Molecular 
mechanisms controlling mutation rates are 
themselves indirectly subject to natural selec- 
tion through genetic modifications they pro- 
duce (second-order selection) (/, 2). The 
linkage between selected mutations and the 
alleles responsible for their generation is par- 
ticularly high in bacteria because their gene- 
transfer and recombination rates are generally 
low. Consequently, when adaptation is limit- 
ed by the supply of mutations, selection was 
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shown to favor strains having constitutively 
increased mutation rates. Such strains display 
high mutation rates owing to the loss of 
genetic fidelity functions, e.g., mutational in- 
activation of the mismatch repair system re- 
sults in a 10?- to 10*-fold increase in mu- 
tagenesis (3). The selection of constitutive 
mutators and their role in adaptive evolution 
of bacteria has been supported by in vivo and 
in vitro experimental evolution, computer 
modeling, molecular evolution, and studies of 
natural isolates (/, 2). 

Mutation rates in bacteria can also be in- 
creased by stress-induced reversible activation of 
some gene functions, which results in a transient 
mutator phenotype, the SOS response being a 
paradigm of such a process (4). However, the 
evolutionary significance of stress-inducible mu- 
tagenesis in bacterial evolution remains a subject 
of intense debate (5, 6). While some argue that it 
is a consequence of a genetically programmed 
evolutionary strategy which, by increasing mu- 
tagenesis, increases the probability of generation 
of adaptive variants, others argue that mutations 
arise in stressed bacteria only as an accidental 
consequence of accumulation and/or processing 
of DNA lesions. However, these hypotheses are 
based on results obtained with laboratory strains. 
It is difficult to assess the evolutionary signifi- 
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cance of any phenomenon without knowing its 
frequency and ecological distribution in natural 
populations, as well as their physiological and 
genetic determinants. With this premise, we 
have studied stress-induced mutagenesis pheno- 
types among 787 worldwide natural isolates of 
Escherichia coli from diverse ecological niches: 
commensal and pathogenic isolates from a vari- 
ety of hosts and isolates from air, water, and 
sediments (7). To mimic stress conditions com- 
monly encountered by bacteria in natural envi- 
ronments, we used progressive starvation fol- 
lowing an exponential growth phase occurring in 
colonies. We chose colonies, instead of liquid 
cultures, because the primary natural FE. coli 
habitat is the gut of warm-blooded animals, 
where it can be found in tightly packed commu- 
nities. In secondary environments, like soil and 
water, bacterial cells also tend to aggregate and 
form (micro)colonies and biofilms. 

Diversity of constitutive and colony- 
aging induced mutation rates among nat- 
ural isolates of EF. coli. To estimate mu- 
tagenesis in aging colonies (MAC) of natural 
isolates of E. coli, we measured the frequency 
of mutations conferring resistance to rifam- 
picin (Rif®) in 1-day- (D1) and 7-day- (D7) 
old colonies (7). For all strains, the median 
values of the frequency of Rif® mutations 
were 5.8 X 10°? on D1 and 4.03 X 107 on 
D7 (Fig. 1, A and B). Thus, the frequency of 
mutations increased on average sevenfold be- 
tween D1 and D7 [(Fig. 1C) Mann-Whitney: 
P < 0.0001], while the median number of 
colony-forming units (CFU) increased 1.2- 
fold. The E. coli K12 MG1655 laboratory 
strain showed a 5.5-fold increase in frequen- 
cy of Rif® mutagenesis and a 1.7-fold in- 
crease in CFU. The increase in CFU from D1 
to D7 was not correlated with the increase in 
the D7/D1 ratio of mutation frequency. 

Strains having D1 mutation frequencies 
>10-fold or >100-fold higher than the median 
D1 mutation frequency of all the strains repre- 
sented 3.3 and 1.4% of isolates, respectively, 
which corresponds to previous reports on the 
frequency of constitutive mutators in natural E. 
coli populations (8—/0). The D7/D1 mutation 
frequency ratio showed that 40% of strains had 
more than a 10-fold, and 13% more than a 
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100-fold, increase in mutagenesis over 7 days 
(Fig. 1C). Interestingly, constitutive mutagenesis 
and MAC showed a negative correlation [Fig. 
1D (Spearman’s correlation —0.347, P < 
0.0001)]. Hence, strains with a lower DI mu- 
tagenesis show significantly higher MAC than 
those with a higher D1 mutagenesis and vice 
versa (Fig. 1E). This was true even when strong 
constitutive mutators (>10-fold the median val- 
ue of D1 mutagenesis for all strains) were ex- 
cluded from the analysis (//). 

Because resistance to rifampicin is conferred 
by mutations in only one gene, rpoB, we have 
also measured mutagenesis at several other 
genomic targets by using /ac/ papillation assay 
and by scoring mutations conferring resistance 
to other antibiotics in a subset of studied strains 
[see (7) and supporting online text]. The results 
showed that mutagenesis increased genome- 
wide, in a large fraction of natural isolates, as a 
consequence of the stresses they encounter in 
aging colonies. 

Bacterial phylogeny, ecological diver- 
sity, and variability of MAC phenotypes. 
The response of bacteria to stress in aging col- 
onies may depend on the frequency of such 
stresses encountered in the course of their 
evolutionary histories. Because no significant 
correlation was observed between MAC and 
phylogenetic groups of the 341 strains with 
established phylogenies (//), and because phy- 
logeny reflects long-term evolutionary histories, 
we concluded that MAC variability reflects 
recent events. Strains belonging to the same 
phylogenetic group are found often in different 
ecological niches and, thus, face different selec- 
tive pressures. Furthermore, strains belonging to 
different phylogenetic groups can share the 
same ecological niches, for example, the human 
commensals used in our study (/2). 

To determine whether MAC variability re- 
flects differences in the ecological niche, we 
measured the MAC for commensal strains iso- 
lated from different hosts with different diets. 
The results indicated that the host’s nutrition 
might explain some of the variation of MAC 
(Fig. 2, A and B). For example, strains isolated 
from different classes of carnivores [mammals 
(n = 47), birds (n = 7), and reptiles (n = 9)] do 
not show any significant difference between 
their MACs. However, these strains are signifi- 
cantly different from most strains isolated from 
noncarnivorous hosts within the same host class. 
Furthermore, strains isolated from humans do 
not have significantly different MAC from iso- 
lates found in other omnivores. 

We also compared MAC phenotypes of com- 
mensal strains isolated from feces of healthy 
volunteers and pathogenic isolates from patients 
with two main types of pathologies caused by E. 
coli, diarrhea and urinary tract infections (Fig. 2, 
C and D). The results show that both groups of 
pathogens from our collection have similar D7/ 
D1 mutation ratios and these are significantly 
lower than those observed for commensal strains 
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(Mann- Whitney: P < 0.0001). The pathogens 
not only showed significantly lower MAC, but 
also had significantly higher constitutive muta- 
tion rates than commensals (Mann-Whitney: 
P = (0.026). Hence, the causes and consequences 
of the MAC phenotype may be relevant for our 
understanding of infectious diseases. 

The differences in recent evolutionary se- 
lective pressures experienced by bacteria in 
different habitats may be responsible for the 
variability of MAC phenotypes. Other pa- 
rameters (e.g., migration, population size of 
the bacteria or of the host, and composition of 
gut microflora), which were not taken in ac- 
count, may also play an important role in 
shaping MAC phenotypes. 

Physiological requirements for MAC. 
Ten randomly chosen isolates with strong 
MAC phenotype were cultivated under differ- 
ent experimental conditions, and the mutagen- 
esis generating Rif® phenotype was measured 
(Table 1). We found that only enrichment of 
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869 plates with a carbon source, e.g., galactose 
(or arabinose; not shown) significantly reduced 
MAC. The addition of amino acids, nitrogen, or 
phosphate sources had no effect. No signif- 
icant increase in the D7/D1 mutagenesis 
ratio was observed when bacteria were kept 
in liquid 869 medium. Therefore, carbon- 
source starvation is necessary but not suf- 
ficient for the increase in MAC. 

An important difference between structured 
(colony) and nonstructured (liquid culture) envi- 
ronments is changing gradients of various chem- 
ical and physical parameters. Bacterial colonies 
are nearly completely anaerobic during the first 
days of growth because of intense metabolic 
activity of the bacteria on the surface of the 
colony. When nutrients are exhausted, metabolic 
activity slows, and oxygen penetrates the colony 
(/3). Under anaerobic conditions, the D7/D1 
ratio of mutagenesis to Rif® was only 1.9 (Table 
1), suggesting that oxidative metabolism is re- 
quired for MAC. Other works suggest that 
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Fig. 1. Polymorphism of constitutive mutagen- 
esis (D1) and mutagenesis in aging colonies (D7) 
among E. coli natural isolates. Histograms (A), 
(B), and (C) represent the distribution of log- 
transformed D1, D7, and D7/D1 mutator fre- 
quencies, respectively. Mutator phenotypes 
have been estimated by measuring the capacity 
of different strains to generate mutations con- 
ferring rifampicin resistance. (D) The negative 


correlation between log-transformed D7/D1 ratios and log-transformed D1 mutator phenotype 
(Spearman's correlation —0.347, P < 0.0001). Each open circle represents one strain. (E) Median 
values of D1 (filled gray circles) and of D7 (filled black circles) mutagenesis are represented for each 
strain. The strains were ordered according to ascending D7/D1 mutation ratios for each strain. This 
illustrates the origin of the negative correlation shown in (D), and reveals an apparent mutual 
exclusion of high constitutive and high MAC mutation phenotypes. 
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“adaptive” mutagenesis in starving E. coli K12 is 
suppressed under anaerobic conditions, or by 
addition of oxygen scavengers, but increased in a 
superoxide dismutase mutant (/4). Together 
these observations and our data suggest that 
oxidative stress is a major contributor to MAC. 

Genetics of MAC: C4750, a case study. 
To study the molecular mechanisms of the 
MAC phenomenon, we chose one natural iso- 
late with a strong MAC phenotype as a model 
(Tables 2 and 3). Although one strain cannot 
represent the genetic diversity of the species, 


this strain responded just as the other nine 


5 1.8 ’ strains did when it was exposed to different 
eS ZG 6) d x a AO oO experimental conditions described in Table 1. 
oe 14 Q Fo Mutational spectra for MAC of C4750 were 
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Fig. 2. Ecological diversity and the variability of mutator phenotypes among £. coli natural isolates. 
(A and B) Fecal isolates were classified according to their host diet; carnivora- (n = 63), insectivora- 
(n = 39), granivora- (n = 33), hematovora- (n = 4), herbivora- (n = 38), omnivora- (n = 32), 
human commensals- (n = 217), nectarivora-, (n = 47), piscivora- (n = 21), frutivora- (n = 9), and 
zooplankton-eating animals (n = 3). (C and D) The role of bacterial life-style. Human strains are 
classified according to their commensal or pathogenic origin: commensal (n = 217), enteroinvasive 
[Shigella (n = 73), and enteroinvasive E. coli (n = 13)], and urinary tract infections (UTI; n = 91) 
strains. (A) and (C) represent log-transformed constitutive mutagenesis (D1; filled gray circles) and 
mutagenesis in aging colonies (D7; filled black circles), whereas (B) and (D) show log-transformed 
D7/D1 mutagenesis ratio (open circles). Mutator phenotypes have been estimated by measuring 
the capacity of different strains to generate mutations conferring rifampicin resistance. Data points 


represent the pooled mean values (+SE) for a given group of isolates. 


Table 1. The physiological requirements of mutagenesis in aging col- 
onies. Mutagenesis to rifampicin resistance in aging colonies of 10 
randomly chosen bacterial strains was measured under different experi- 
mental conditions. Pooled values of Day 1 and Day 7 mutation fre- 
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e spectively) were G:C—>A:T transitions. The in- 
crease in the frequency of these mutations be- 
tween D1 and D7 is significant (x* test: P = 
0.03). In addition, we sequenced rpoB genes 
from 90 Rif® mutants belonging to 23 different 
natural isolates having a strong MAC and found 
that, similar to C4750, G:C—A:T transitions 
represent the majority of substitutions (57 and 
82% on D1 and D7, respectively; x? test: P = 
0.009). However, unlike C4750, those strains 
also showed a significant increase of G:C>T:A 
transversions (0 and 11% on D1 and D7, respec- 
tively; x7: P = 0.018). 


quency obtained for each exprimental condition for 10 strains are 
presented. 7 Mutation frequency significantly higher relative to mutation 
frequency on 869 plate. | Mutation frequency significantly lower relative 
to mutation frequency on 869 plate. 
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ean median value Fa median value . 
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(238) felatve tO (relative to 869 feline te (relative to 869 
869 plate 869 plate 
plate) plate) 
869 plate 3.0 + 1X 10-8 i 3.2 + 0.7 X 10-6 A - 
Anaerobic 2.2 + 08x 1077 71 <0.0001 ¢ 4.3+2.2x 10-7 0.08 <0.0001 | 
869 plate 
869 plate 3.6 + 13x 1078 1.2 0.46 9.2 +3.3 x 108 0.02 <0.0001 | 
1% galactose 
869 plate 9.8 + 1.7 107° 0.3 0.42 44+ 1.4 x 107° 2.1 0.19 
1% casamino 
acids 
869 plate Vf 110-8 0.5 0.36 1.44 0.3 X 10-6 0.6 0.52 
100 mM NH,ClL 
869 plate 16 + 0.2 x 107° 0.5 0.29 3.2+0.5 x 10° 1.8 0.13 
10 mM K,HPO, 
869 liquid 5.8 + 2.1x 108 3.7 <0.0001 7 2.0 + 0.6 X 1077 0.03 <0.0001 | 
*Significance threshold: P = 0.05. 
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Because the G:C—A:T transitions repre- 
sent the majority of mutations on D1 and D7 in 
the C4750 strain, we used competition experi- 
ments to test whether rpoB Rif® alleles carrying 
these mutations confer a growth advantage in 
old colonies, which would bias the observed 
mutation frequency. The results suggested that 
there is no growth advantage and that the in- 
crease in the frequency of G:C—A:T mutations 
in rpoB is caused by increased mutagenesis (see 
supporting online text). 

(ii) We measured mutagenesis in lacZ al- 
leles in colonies incubated 7 days on agarose 
plates (7). On agarose, C4750 showed a 3.8- 
fold increase in frequency of Rif® mutagenesis 
(Mann-Whitney: P < 0001) and a 3-fold 
decrease in CFU. Four reversions in the lacZ 
gene showed a significant increase over 7 
days: G:C A:T transitions (9.9-fold; Mann- 
Whitney: P = 0.04) and —1G (10-fold; 
Mann-Whitney: P = 0.01), —2CG (681-fold; 


Table 2. Genetic analysis of mutagenesis in aging colonies of the E. coli 
C4750 natural isolate. For the change of number of viable cells in aging 
colonies of different C4750 derivatives see supporting online text. | Mutant 


Mann-Whitney: P = 0.02), and —1A frame- 
shifts (97-fold; Mann-Whitney: P = 0.02). 

The most frequent MAC substitution mu- 
tations in rpoB and lacZ genes are G:C—> A:T 
transitions that are increased in many E. coli 
DNA-repair mutants [e.g., ung”, nth”, nei, 
ada~ and ogt’, mutS~ (3)], so it was not 
possible to assign a specific mutagenic path- 
way for MAC. However, a coincidently 
high increase of frameshift mutations in the 
lacZ gene suggests the involvement of a 
mismatch-repair deficiency or insufficien- 
cy (/5) [see (iv) below]. 

(iti) Because in C4750 strain, MAC is asso- 
ciated with carbon-source starvation and oxida- 
tive shock, we investigated different pathways 
contributing to bacterial survival under such 
stresses. The o* subunit of RNA polymerase is 
the regulator of expression of genes involved in 
responses to diverse stresses including station- 
ary phase, starvation, osmotic, acid, heat, and 
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oxidative shocks (/6). The inactivation of ei- 
ther of rpoS or hfg, which encode transcription 
factor o* and the protein controlling rpoS 
mRNA translation, respectively (/6), abolished 
the increase in MAC on D7 (Table 2). The 
inactivation of rssB, resulting in high constitu- 
tive levels of o* (J6), led to significantly in- 
creased MAC over the parental strain level, 
further confirming positive control of MAC by 
the RpoS regulon (Table 2). 

We examined the role of the cyaA and crp 
genes (encoding adenylate cyclase and cAMP 
receptor protein, respectively) in MAC be- 
cause MAC can be attenuated by adding 
more sugar to the plates and because carbon- 
source starvation results in activation of the 
cAMP/CRP regulatory network (17). It ap- 
pears that MAC requires the activity of both 
genes (Table 2). 

(iv) The spectra of mutations in the rpoB 
and /acZ genes implies the involvement of the 


strain mutation frequency significantly higher relative to mutation frequency 
of parental strain. | Mutant strain mutation frequency significantly lower 
relative to mutation frequency of parental strain. 
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Frequency of mutagenesis to Rif® on day 7 


Genshine Ratio of P value of Ratio of P value of 
yP Mean median values Mann-Whitney test* Mean (-SE) median values Mann-Whitney test* 
(SE) mutant/ Mutant versus ~ mutant/ Mutant versus 
parental strain parental strain parental strain parental strain 
Parental 3.5 + 0.9X 10-8 1 = 2.7 + 0.6 X 10°© 1 = 
strain 
rpoS 8.6 + 4x 10-9 0.2 0.028 | 32.2 10-8 0.005 <0.0001 | 
hfq 1.2 + 06x 1078 0.4 0.16 9.3 + 6x108 0.004 0.0007 | 
rssB 13 + 04x 1078 0.7 0.24 6.8 + 1x 10-6 5.5 0.0012 
cyaA 1.2 + 05x 10-8 0.4 0.10 3.4 + 2x 108 0.009 <0.0001 | 
crp 25 + 1x108 0.7 0.71 7.7 4+ 1x 1078 0.06 0.007 | 
recA 16 + 04x 1078 0.6 0.29 2.1 + 0.6 X 1077 0.05 <0.0001 | 
lexA1 2.4 + 0.6 x 10-8 1.4 0.47 18 + 0.5 x 10° 0.9 0.82 
polB 3.4 + 09x 10-8 1.3 0.28 6.0 + 2x10’ 0.3 0.0007 | 
uvrA 7.2 + 3x 108 2.9 0.12 2.9 + 1X 10-& 0.6 0.98 
polA 1.5 + 08x 1077 48 0.045 {+ 25 + 2x10 4.1 0.12 
recB 2.3 + 0.7 Xx 10-8 0.9 0.99 3.4 + 1X 10° 1.1 0.75 
mutS 3.7 + 0.7X 107° 195 <0.0001 T 18 + 03x 10~° 1.25 0.96 


*For this nonparametric test, all values obtained from at least three independent experiments (each with three independent cultures) for each genotype were used (significance 
threshold, P = 0.05). 


Table 3. The role of MutL and MutS protein overproduction in mutagenesis 
in aging colonies of the F. coli C4750 natural isolate The strains designated 
pvector, pmutL*, and pmutS* were C4750 derivaties carrying pACYC184, 
pMQ341, and pMQ339 plasmids, respectively. pvector was considered the 


parental strain and was used for comparison with the other two strains carrying 
plasmids. { Mutant strain mutation frequency significantly higher relative to 
mutation frequency of parental strain. | Mutant strain mutation frequency 
significantly lower relative to mutation frequency of parental strain. 


Frequency of mutagenesis to Rif® on day 1 Frequency of mutagenesis to Rif® on day 7 


P value of P value of 


Genotype Ratio of median Ratio of median 


Mean values imutant/ Mann-Whitney test* Mean values trutant/ Mann-Whitney test* 
(SE) : Mutant versus (+SE) : Mutant versus 
parental strain ‘ parental strain : 
parental strain parental strain 
Pvector 9.0 + 3x10? 1 - 13 4 1} - 
0.3 x 1077 
PmutL* 9.4 42x10 1.39 0.84 85 = 0.70 0.099 
0.8 x 107-8 
PmutS* if 0.85 0.81 2:2. 0.17 0.006 | 
0.8 X 10°? 0.5 x 10-8 


*For this nonparametric test, all values obtained from at least three independent experiments (each with three independent cultures) for each genotype were used (significance 
threshold, P = 0.05). 
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mismatch-repair system in MAC. The mis- 
match-repair system is a major contributor to 
DNA replication and recombination fidelity in 
growing E. coli K12 cells, but seems to be 
nonfunctional during stationary phase (/8). 
During stationary phase the mutS gene tran- 
script and the MutS [mismatch recognizing and 
binding (4)] protein decrease to undetectable 
levels through a RpoS- and Hfq-dependent 
mechanism (/9). Because MAC in C4750 
strain was also found to be rpoS and hfq 
dependent, we looked for the potential involve- 
ment of mismatch-repair down-regulation in 
MAC. The mutS mutant did not show any 
significant elevation of MAC (Table 2). By 
contrast, overproduction of the MutS protein 
significantly depressed MAC, whereas over- 
production of the MutL protein had no effect 
(Table 3). Hence, down-regulation of mismatch 
repair seems to be responsible for the increase 
in MAC, and consequently, the majority of 
these mutations arise during DNA synthesis. 

(v) It has been previously reported that E. 
coli K12 MAC and “adaptive” mutagenesis de- 
pend on the induction of the SOS system (20, 
21). The SOS regulon is composed of at least 40 
genes (most of which encode DNA-repair func- 
tions) and is induced by a variety of stresses that 
damage DNA and/or block DNA replication 
(22). The RecA protein acts as a positive regu- 
lator, whereas LexA is the repressor of the SOS 
genes (4). Our data indicate that 95% of muta- 
tions that arise in aging colonies were RecA 
dependent and LexA independent (Table 2). We 
tested two other SOS genes, wvrA and polB, 
encoding the nucleotide excision and DNA poly- 
merase II (Pol II) proteins, respectively. The 
inactivation of the wvrA gene had no significant 
effect, whereas 73% of mutations in aging col- 
onies were polB dependent (Table 2). The error- 
prone activity of this polymerase has also been 
observed in other studies (23, 24). Because Pol 
II has been shown to be involved in the response 
to oxidative damage (25), and because MAC is 
oxygen dependent (Table 1), it may be that this 
polymerase participates in the error-prone copy- 
ing of oxidative lesions. 

The inactivation of the po/A gene (not a part 
of the SOS regulon), which encodes DNA poly- 
merase I [involved in gap-filling steps during 
DNA replication and nucleotide excision repair 
(4)], did not modify the MAC phenotype (Table 
2). The involvement in MAC of at least two 
SOS genes, recA and polB, suggests that either 
the basal levels of expression of these genes is 
sufficient for MAC or that some SOS genes may 
be regulated in a noncanonical (LexA-indepen- 
dent) manner in quiescent cells. SOS induction 
and “adaptive” mutagenesis in E. coli K12 qui- 
escent cells has been shown to depend on the 
recA and recB genes (26, 27), suggesting the 
involvement of double-strand breaks and re- 
combination in that mutagenic pathway. How- 
ever, recB gene function is not required for 
MAC of the C4750 strain (Table 2). 
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Fig. 3. Probability of fixation of mutator alleles as a function of their effect on mutation rate. The 
adaptation of a population of 10° cells was simulated in environments in which stress phases occurred 
(on average) every 33 generations and lasted for (on average) 10 generations. In this environment, the 
overall fraction of time spent under stress (K) is 23%. (A) Constitutive mutator alleles (circles) had a 
permanent X-fold increase in mutation rate, whereas the stress-inducible mutators (triangles) had an 
X-fold increase in mutation rate only during stress phases. (B) The strength of the stress-inducible 
mutators has been normalized to reflect its mean effect on the generation of mutations over stress and 
nonstress periods [X’ = K X X+(1 — K) X 1]. Each point is the fixation probability over 500 simulations. 
Error bars represent the standard error of the fixation frequency over 100 simulations. 


Simulation of the selection of stress-in- 
duced mutator alleles and their effects. We 
used computer simulations previously developed 
to study the role of constitutive mutators in the 
adaptive evolution of bacteria (7, 28). We as- 
sumed that stress-inducible mutator alleles can- 
not be selected for or against because of their 
pleiotropic effects. We did not take into account 
mutations that result in immediately improved 
survival to these stresses, but only those that 
have an impact on the longer-term adaptation of 
bacterial populations. The simulations showed 
that stress-induced mutator alleles can be select- 
ed almost as efficiently as constitutive mutators 
[(Fig. 3) and supporting online text], which is in 
contrast to the observed higher proportion of 
MAC mutators relative to constitutive mutators 
among natural isolates (for 10-fold mutators, 40 
versus 3.3%, respectively). However, the wide 
distribution of MAC phenotypes among natural 
isolates (Fig. 1) and a variety of genetic mecha- 
nisms involved (/4) suggest that the number of 
genes involved in the MAC phenotype is higher 
than the number of genes involved in the consti- 
tutive mutator phenotype. Consequently, the se- 
lection of MAC mutator genotypes in nature 
could be more efficient than the selection of 
constitutive ones under directional selection. 

The selection of stress-inducible mutators in 
nature may be also due to pleiotropic effects of 
strongly selected functions, e.g., DNA repair or 
other survival-promoting systems. However, re- 
gardless of the nature of selection of stress- 
inducible alleles, mutations produced under 
stress could represent a large proportion of over- 
all mutations and may have evolutionary conse- 
quences. Since the mean strength of MAC 
among natural isolates is about 10-fold, if phases 
of stress were to occur 25% of the time, the 


proportion of mutations produced under stress 
would be, on average, 46%. Because some of 
those mutations might be adaptive, stress- 
inducible mutagenesis would increase the speed 
of adaptation. Our simulations showed that a 
population with a 10-fold stress-inducible muta- 
tor phenotype will adapt up to 15% faster than a 
nonmutator population, whereas a population 
having a 100-fold stress-inducible mutator phe- 
notype will adapt up to 38% faster. 

Increase in adaptation rate due to stress-in- 
ducible mutagenesis may also limit the selection 
of constitutive mutator alleles. The fixation of an 
allele improving the rate of adaptation has been 
shown to decrease the selection for other alleles 
improving the rate of adaptation (29). Indeed, 
our simulations showed that the frequency of 
fixation of a constitutive mutator decreases as a 
function of increased stress-inducible mutagene- 
sis (Fig. 4). The relative decrease in the fixation 
of a constitutive mutator allele was almost per- 
fectly correlated with the relative improvement 
of the adaptation rate due to stress-inducible 
mutagenesis. This effect of MAC on adaptation 
rate could explain the negative correlation ob- 
served between constitutive mutagenesis and 
MAC (Fig. 1D). 

Concluding remarks. Most natural isolates 
of E. coli exhibits increased mutation rates under 
stress encountered in aging colonies (Fig. 1). 
MAC is characteristic to each strain, varies 
greatly, and increases as a consequence of stress 
imposed by carbon-source starvation and oxida- 
tive shock (Table 1). Because this stress-induc- 
ible mutagenesis is genetically controlled (Table 
2), and therefore subject to selective pressures, 
its high variability should reflect the diversity of 
selective pressures in different environments. In- 
deed, our simulations showed that the strength of 
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Fig. 4. Probability of 
fixation of a 100-fold 
constitutive mutator 
in populations with 
various strengths of 
stress-inducible | mu- 
tagenesis. Populations 
evolved under the 
conditions described 
in Fig. 3. During stress 
periods the nonmuta- 
tor population had an 
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inducible mutagene- 
sis), whereas the 100- 
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tator cells had the strongest increase in mutation rate, between 100- and X-fold. During nonstress 
periods only constitutive mutator individuals had a 100-fold increased mutation rate. Each point is 
the fixation probability over 1000 simulations; error bars represent the standard error of the 


fixation frequency over 100 simulations. 


the selected stress-induced mutator alleles is pos- 
itively correlated with the strength of selection 
and negatively with the frequency of such stress- 
es (11). The fact that the ecological niche from 
which the strain was isolated appears to be the 
major determinant of mutator phenotype (Fig. 
2), and the absence of correlation between bac- 
terial phylogeny and mutator phenotype suggest 
that the evolution of stress-inducible mutation 
rates is rapid on an evolutionary time scale. This 
hypothesis is corroborated by our simulations 
where the selection of stress-induced mutator 
alleles occurred on a range of hundreds to thou- 
sands of generations. 

What are the mechanisms responsible for 
these variable mutator phenotypes? Computer 
simulations, which took into account only the 
mutagenic effect, clearly show that second- 
order selection may drive the evolution of stress- 
inducible mutator alleles (Fig. 3). However, our 
results also agree with the possibility that stress- 
induced mutagenesis is a side effect of some 
stress-resistance mechanisms emerging from 
first-order selection. For example, we have ob- 
served that the down-regulation of mismatch 
repair seems to be involved in the increase of 
mutation rates in aging colonies (Table 3). Tran- 
sient limitation of mismatch-repair activity was 
also reported to be involved in “adaptive” mu- 
tagenesis (30). Because such down-regulation 
was shown to be RpoS dependent (19), and 
because the RpoS-regulated network normally 
increases stress resistance (/6), it may be that 
this down-regulation was selected for to avoid 
the energy cost incurred by the repair functions 
when food is limited. Moreover, potentially mu- 
tagenic SOS polymerases are known to increase 
survival of bacteria exposed to DNA-damaging 
agents, while increasing mutagenesis (3/, 32). In 
aging colonies, DNA lesions may be processed 
by such error-prone repair mechanisms, as sug- 
gested by the implication of Pol II polymerase in 
MAC (Table 2) and Pol IV in “adaptive” mu- 
tagenesis (33). In both cases, the resulting pro- 
duction of deleterious mutations may have a less 
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severe immediate effect on fitness than the direct 
cost imposed by efficient error-free survival 
mechanisms. It is plausible that different stress- 
resistance mechanisms have different pleiotropic 
effects, and thus different mutator phenotypes 
may be selected in different ecological niches. 
However, independent of the selective pressure 
that led to the emergence of MAC, genetic di- 
versity produced during stress seems to be a 
substrate for natural selection. Consequently, 
MAC is expected to have an important impact 
on bacterial adaptive evolution. 

It has been proposed that natural bacterial 
populations switch between constitutive mutator 
and nonmutator genotypes by mutations inacti- 
vating mismatch-repair genes and by horizontal 
gene transfer that restores their function (34). 
This study shows that bacteria switch phenotyp- 
ically between high and low mutation rates de- 
pending on environmental conditions. This phe- 
notypic switch seems to depend on a variety of 
molecular mechanisms, as suggested by the vari- 
ability of MAC phenotypes among natural iso- 
lates and different genetic requirements for 
MAC in different E. coli strains (14, 26, 35). 
Because MAC-related mechanisms have been 
observed both in bacteria [e.g., Mycobacterium 
smegmatis (36) and Pseudomonas putida (37)]| 
and yeast (38), we are aware of the possibility 
that the age-related increase in cancer incidence 
(39), somatic mutations (39), and germline mu- 
tations (40) could also be due to the accumula- 
tion of mutations in quiescent stem cells by 
mechanisms akin to MAC. 


References and Notes 

. A. Giraud, M. Radman, |. Matic, F. Taddei, Curr. Opin. 
Microbiol. 4, 582 (2001). 

2. O. Tenaillon, F. Taddei, M. Radman, |. Matic, Res. 
Microbiol. 152, 11 (2001). 

. J. H. Miller, Mutat. Res. 409, 99 (1998). 

4. E. C. Friedberg, G. C. Walker, W. Siede, DNA Repair 
and Mutagenesis (American Society for Microbiology, 
Washington, DC, 1995). 

5. R. E. Lenski, P. D. Sniegowski, Science 269, 285 
(1995). 

6. M. Chicurel, Science 292, 1824 (2001). 


es 


WwW 


RESEARCH ARTICLES 


7. Materials and methods are available as supporting 
material on Science Online. 

8. J. E. LeClerc, B. Li, W. L. Payne, T. A. Cebula, Science 
274, 1208 (1996). 

9. |. Matic et al., Science 277, 1833 (1997). 

10. E. Denamur et al., J. Bacteriol. 184, 605 (2002). 

11. |. Bjedov et al., data not shown. 

12. P. Duriez et al., Microbiology 147, 1671 (2001). 

13. A. Peters, J. Wimpenny, J. Coombs, /. Gen. Microbiol. 
133, 1257 (1987). 

14. B. A. Bridges, Mutat. Res. 408, 1 (1998). 

15. C. G. Cupples, M. Cabrera, C. Cruz, J. H. Miller, Ge- 
netics 125, 275 (1990). 

16. R. Hengge-Aronis, Microbiol. Mol. Biol. Rev. 66, 373 
(2002). 

17. M.S. Saier-Jr, T. M. Ramseier, J. Reizer, in Escherichia coli 
and Salmonella Cellular and Molecular Biology, F. C. 
Neidhardt et al., Eds. (American Society for Microbiolo- 
gy, Washington, DC, 1996), vol. 1, pp. 1325-1343. 

18. D. Brégeon, |. Matic, M. Radman, F. Taddei, J. Genet. 
78, 21 (1999). 

19. H.-C. T. Tsui, G. Feng, M. E. Winkler, J. Bacteriol. 179, 
7476 (1997). 

20. G. J. McKenzie, R. S. Harris, P. L. Lee, S. M. Rosenberg, 
Proc. Natl. Acad. Sci. U.S.A. 97, 6646 (2000). 

21. F. Taddei, |. Matic, M. Radman, Proc. Natl. Acad. Sci. 
U.S.A. 92, 11736 (1995). 

22. J. Courcelle, A. Khodursky, B. Peter, P. O. Brown, P. C. 
Hanawalt, Genetics 158, 41 (2001). 

23. |. Tessman, M. A. Kennedy, Genetics 136, 439 (1994). 

24. R. Napolitano, R. Janel-Bintz, J. Wagner, R. P. Fuchs, 
EMBO J. 19, 6259 (2000). 

25. M. Escarceller et al., /. Bacteriol. 176, 6221 (1994). 

26. F. Taddei, J. A. Halliday, M. Matic, M. Radman, Mol. 
Gen. Genet. 256, 277 (1997). 

27. R. S. Harris, S. Longerich, S. M. Rosenberg, Science 
264, 258 (1994). 

28. O. Tenaillon, B. Toupance, H. Le Nagard, F. Taddei, B. 
Godelle, Genetics 152, 485 (1999). 

29. O. Tenaillon, H. Le Nagard, B. Godelle, F. Taddei, Proc. 
Natl. Acad. Sci. U.S.A. 97, 10465 (2000). 

30. R. S. Harris et al., Genes Dev. 11, 2426 (1997). 

31. S. G. Sedgwick, P. A. Goodwin, Proc. Natl. Acad. Sci. 
U.S.A. 82, 4172 (1985). 

32. L. J. Runyen-Janecky, M. Hong, S. M. Payne, Infect. 
Immunol. 67, 1415 (1999). 

33. G. J. McKenzie, P. L. Lee, M. J. Lombardo, P. J. Hastings, 
S. M. Rosenberg, Mol. Cell 7, 571 (2001). 

34. E. Denamur et al., Cell 103, 711 (2000). 

35. P. L. Foster, Annu. Rev. Genet. 33, 57 (1999). 

36. P. Karunakaran, J. Davies, /. Bacteriol. 182, 3331 
(2000). 

37. S. Saumaa, A. Tover, L. Kasak, M. Kivisaar, J. Bacteriol. 
184, 6957 (2002). 

38. D. F. Steele, S. Jinks-Robertson, Genetics 132, 9 
(1992). 

39. R. Holliday, Understanding Ageing (Cambridge Univ. 
Press, Cambridge, 1995). 

40. S. S. Sommer, W. A. Scaringe, K. A. Hill, Mutat. Res. 
487, 1 (2001). 

41. We thank C. Amorin, J. P. Coutanceau, P. Duriez, L. 
Espinoza, E. Mancera, M. Selva, and A. Valera for tech- 
nical help; |. Le Gall for DNA sequencing; B. Picard for 
providing isolates from humans; and A. Giraud, L. Le 
Chat, A. Lindner, G. Simhen, and E. Stewart for critical 
reading of the manuscript. This work was supported by 
grants from the Programme de Recherche Fondamen- 
tale en Microbiologie et Maladies Infectieuses et Para- 
sitaires-MENRT and “Programme Environnement et 
Santé-MATE.” The establishment of the collection of E. 
coli isolates from animal hosts was supported by grants 
from CONACyT 0028 and DGAPA IN208601. |.B. was 
supported by a fellowship from the Association pour la 
Recherche sur le Cancer. 


Supporting Online Material 
www.sciencemag.org/cgi/content/full/300/5624/1404/ 
DC1 

Material and Methods 

SOM Text 

References 


10 January 2003; accepted 4 April 2003 


1409 


